Neutron diffraction evidence of microscopic charge inhomogeneities in the CuC>2 plane 

of superconducting La 2 - x Sr x Cu04 (0 < x < 0.30) 
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We present local structural evidence supporting the presence of charge inhomogeneities in the CuC>2 
planes of underdoped La2-^Sr a; Cu04. High-resolution atomic pair distribution functions have been 
obtained from neutron powder diffraction data over the range of doping < x < 0.30 at 10 K. 
Despite the average structure getting less orthorhombic we see a broadening of the in-plane Cu-0 
bond distribution as a function of doping up to optimal doping. Thereafter the peak abruptly 
sharpens. The peak broadening can be well explained by a local microscopic coexistence of doped 
and undoped material. This suggests a crossover from a charge inhomogeneous state at and below 
optimal doping to a homogeneous charge state above optimal doping. The strong response of the 
local structure to the charge-state implies a strong electron-lattice coupling in these materials. 
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There is mounting interest in the possibility that the 
charge distribution in the CuC>2 planes of the high- 
temperature superconductors is microscopically inhomo- 
geneous |l],|2j and that this has a bearing on the high 
temperature superconductivity itself. An inhomogeneous 
charge distribution has been predicted theoretically |j] Q . 
There is also mounting experimental evidence that micro- 
scopic charge inhomogeneities exist in particular cuprate 
samples. The most compelling evidence for this comes 
from the observation from neutron diffraction of stripes 
of localized charge in L^-z-yNdySr^CuO^ |§|§. How- 
ever, such direct evidence for charge stripes has only been 
seen in insulating compounds or in samples where this 
ordering competes with superconductivity ||. The ev- 
idence supporting the presence of dynamic local charge 
stripe distributions in the bulk of superconducting sam- 
ples is based primarily on the observation of incommen- 
surate spin fluctuations []l0|-^3|. We present local struc- 
tural evidence that supports the fact that the charges are 
inhomogeneous in the underdoped and optimally doped 
region of the La2- a; Sr 2 ;Cu04 phase diagram, consistent 
with the presence of charge domains or dynamic charge- 
stripes. The charge distribution becomes homogeneous 
on crossing into the overdoped region. The clear observa- 
tion of these effects in the local structure also underscores 
the point that there is a strong electron-lattice coupling, 
at least to particular distortion modes, in these materials. 
We note that the disappearance of the structural distor- 
tions, and therefore the charge inhomogeneities, corre- 
lates with the disappearance of the normal-state pseudo- 
gap [Q in these materials. 

The presence of charge inhomogeneities in the CuC>2 
planes implies profound consequences for the local struc- 
ture. It is well known that the average Cu-0 bond length 
changes as the charge state of copper changes. Thus, the 



Cu-0 bond in La2- a; Sr a ;Cu04 shortens from 1.904 A to 
1.882 A as x changes from to 0.2 and the aver- 
age copper charge changes from 2+ to ~ 2.2+. This is 
a generic feature of variably dopable HTS samples and 
comes about because the Cu-0 bond is a covalent anti- 
bonding state which is stabilized by removing electron 
density from it. Clearly, if the doped charge in the Cu02 
planes is inhomogeneously distributed, such that some 
copper sites have more charge than others, a distribution 
of in-plane Cu-0 bond lengths will exist. A high res- 
olution measurement of the in-plane Cu-0 bond length 
distribution as a function of doping will therefore reveal 
the extent of charge inhomogeneities. 

We have used the atomic pair distribution function 
(PDF) analysis of neutron powder diffraction data to 
measure accurate Cu-0 bond length distributions with 
high resolution for a series of La2_ £! ;Sr :E Cu04 samples 
with (0 < x < 0.3). The ability of high-resolution PDF 
studies to reveal local bond-length inhomogeneities which 
are not apparent in the average structure has been clearly 
demonstrated (l6) . The samples studied cover the range 
from undoped, through underdoped and optimally doped 
(x = 0.15) to the overdoped regime. We find that at 
10 K the mean-square width of the bond- length distri- 
bution, a 2 , increases approximately linearly with x until 
optimal doping above which it sharply decreases and re- 
turns to the value of the undoped material by x = 0.25. 
This is strong evidence for charge inhomogeneities in the 
under- and optimally-doped regimes as we discuss below. 
This increase in bond-length distribution can be well ex- 
plained by a linear superposition of the local structures 
of undoped and heavily doped material. 

Samples of Lag-xSr^CuC^ with x = 0.0, 0.05, 0.1, 
0.125, 0.15, 0.2, 0.25, and 0.3 were made using stan- 
dard solid state synthesis. Mixtures of La 2 03, SrC03 
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and CuO were reacted at temperatures between 900°C 
and 1050°C with intermediate grindings, followed by an 
annealing step in air at 1100°C for 100 hours and in 
oxygen at 800°C for 100 hours. The long anneals were 
carried out to ensure doping homogeneity in the sam- 
ples. The lattice c-axis parameter is a sensitive measure 
of oxygen stoichiometry fl2] ]. Its value was determined 
for each of our samples thru Rietveld refinement of the 
data. It varies smoothly with x as expected for stoi- 
chiometric samples and lies on the c-axis vs x curve of 
Radaelli et al. |L5j within the scatter of their data. Neu- 
tron powder diffraction data were collected at 10 K on 
the Special Environment Powder Diffractometer (SEPD) 
at the Intense Pulsed Neutron Source (IPNS) at Argonne 
National Laboratory. Approximately lOg of finely pow- 
dered sample was sealed in a cylindrical vanadium tube 
with He exchange gas. The samples were cooled us- 
ing a closed-cycle He refrigerator. The data were cor- 
rected fl7| for experimental effects and normalized to 
obtain the total structure function S{Q) where Q rep- 
resents neutron momentum transfer. The PDF, G(r) 1 is 
obtained by a Fourier transform of the data according to 
G(r) = | jf °° Q[S(Q) - 1] sin Qr dQ. PDFs from these 
samples are shown in Figs. 5-8 of Ref. [jfjl and in Fig. || 
of this paper. The PDFs examined in this paper used 
data over a range 0.7 A _1 < Q < 28 A -1 . 

We are interested in extracting the width of the dis- 
tribution of in-plane Cu-0 bond- lengths. This informa- 
tion is contained in the width of the first PDF peak at 
~ 1.9 A. The peak-width comes from the thermal and 
zero-point motion of the atoms plus any bond- length dis- 
tribution originating from other effects such as charge 
inhomogeneities. We can determine the latter by consid- 
ering the PDF peak width of this peak as a function of 
doping at 10 K. Three independent measures of the peak 
width all show that the width increases significantly with 
doping up to x = 0.15, beyond which the peak quickly 
sharpens. First, this behavior is evident by simply look- 
ing at the data shown in Fig. |j} This shows the low-r 
region of the PDF around the r = 1.9 A peak as a func- 
tion of x. Since we want to compare the relative peak 
widths (and heights), in this Figure the peaks have been 
shifted to line up the peak centroids and rescaled slightly 
to ensure that the integrated intensity of each peak is 
the same. It is clear that some of the peaks are sig- 
nificantly broader than others with lower peak heights 
and less steeply sloping sides. We have quantified this 
by fitting the peak with a single Gaussian. The Gaus- 
sian is first convoluted with the Fourier transform of the 
step function which was used to terminate the data fll9| . 
This accounts for any termination ripples in the data 
introduced by the Fourier transform of the finite range 
data-set and doesn't introduce any additional parameters 
into the fit. The fitting parameters were peak position, 
scale-factor and width. The baseline is set by the aver- 
age number density of the material and this was fixed at 
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FIG. 1. PDF peak coming from the in-plane Cu-0 bond for 
various doping levels (solid line: undoped case; solid circles: 
overdoped cases; dashed lines: intermediate doping cases). 
The peaks have been shifted so that their centers line up at 
r m = 1.91 A. 

po = 0.07299 A" 3 . The results are shown in Fig. |as the 
solid circles. The mean-square width of the distribution 
increases monotonically (and almost linearly) with x un- 
til x — 0.15. Between 0.15 and 0.2 the peak abruptly 
sharpens and returns to the width of the undoped sam- 
ple by x — 0.25. The same behavior can be obtained 
from the data in a totally model-independent way. If 
the integrated area of a Gaussian is constant the height 
is inversely proportional to the width. Thus, the peak 
height, h, of the rescaled data shown in Fig. [j] should be 
inversely proportional to the width. The open squares in 
Fig. show C/h 2 where h was determined directly from 
the peak maximum in the data and the constant C was 
chosen to make the x = 0.0 points line up. There is ex- 
cellent agreement lending confidence to the results from 
the fitting. 

We would like to discuss possible origins for these dop- 
ing dependent changes in Cu-0 bond-length distribution. 
First we rule out the possibility that it simply comes from 
changes in the orthorhombicity of the sample. The Cu-0 
PDF peak first broadens smoothly with increasing dop- 
ing then dramatically sharpens at a composition close to 
the LTO-HTT structural phase boundary. This behavior 
does not reflect the monotonic decrease in orthorhombic- 
ity of the average structure. Indeed, in the overdoped 



2 




J i I i I i I i I i I i L 



0.05 0.1 0.15 0.2 0.25 0.3 
strontium content, x 

FIG. 2. Peak width of the in-plane Cu-O PDF peak as 
a function of doping obtained by fitting a Gaussian (solid 
circles). The data are plotted as a where a is the Gaussian 
standard deviation. The inverse peak-height-squared of the 
peaks in Fig. 1 scaled so the x = 0.0 points line up is shown 
as open squares. The dashed line is a guide for the eye. 

region the PDF peak returns to the same narrow width 
it had in the undoped material which has the largest 
orthorhombic distortion of any of the samples. We also 
note that the in-plane Cu-0 bonds are not expected to be 
sensitive to the orthorhombic distortion. They lie along 
the unit cell diagonals and not along the unit cell edges 
in the orthorhombic unit cell. In this case an orthorhom- 
bic distortion will change the bond-length but will not 
lead to two distinct bond- lengths. Next we show that 
the observed behavior cannot be explained by doping de- 
pendent changes in the octahedral tilts. The average Jlf| 
(and local) ]l8| tilt amplitude monotonically decreases 
with increasing doping which does not correlate with the 
behavior of the Cu-0 bond length distribution. On the 
other hand, PDF peaks at higher-r do sharpen monoton- 
ically with increasing x reflecting the reduced octahedral 
tilt amplitude |Tsj] . We also rule out the idea that the 
observed Cu-0 peak broadening is an effect of size-effect 
disorder due to doping since the peak sharpens dramat- 
ically above x — 0.2 where the dopant induced disor- 
der should be the greatest. We also note that size-effect 
dopant induced disorder is expected to have a large ef- 
fect on octahedral tilts and a small effect on Cu-0 bond 
lengths since the energy to change a Cu-O-Cu bond angle 



is much less than the energy to stretch the short Cu-0 
covalent bond. Furthermore, the extent of dopant in- 
duced tilt disorder is relatively small in La2- K Sr^Cu04 
as evidenced by the observation that highcr-r PDF peaks 
sharpen on increased doping in this material. These 
peaks sharpen because of a decrease in both the average 
orthorhombicity and octahedral tilt angle with increased 
doping. However, significant size-effect octahedral tilt 
disorder due to the chemical dopants tends to counter 
this effect, as seen in La 2 - x Ba :I ;Cu04 [^0|,^l). We can 
also rule out structural fluctuations associated with the 
HTT-LTO transition. First, these fluctuations will affect 
primarily octahedral tilts and local orthorhombicity (the 
two order parameters of this structural transition) and 
as we have discussed, the Cu-0 bond is expected to be 
quite insensitive to disorder in these parameters. How- 
ever, in addition, we would expect these fluctuations to 
be largest when the structural phase transition temper- 
ature, T s , is closest to our measurement temperature. 
These temperatures are closest for the x = 0.2 sample 
(T s =60 K p2|, T meas =10 K) and this sample exhibits a 
narrow distribution of Cu-0 bond-lengths. The largest 
distribution of Cu-0 bond-lengths is seen for x = 0.15 
where T s = 180 K. 

The observed behavior of the Cu-0 bond-length dis- 
tribution is best explained by the presence of charge in- 
homogeneities. As we have described, the charge-state 
of copper has a direct effect on the Cu-0 bond length 
with the bond-length decreasing with increasing doping. 
Charge inhomogeneities will, thus, give rise to a distri- 
bution of Cu-0 bond lengths. Increased doping will re- 
sult in more Cu-0 bonds being affected and therefore a 
larger measured effect in the PDF, as observed. Above 
optimal doping the PDF peak width abruptly sharpens 
to its value in the undoped material. This is consistent 
with the idea that in the overdoped region, the charge dis- 
tribution in the Cu-0 planes is becoming homogeneous. 

We now discuss independent evidence from the data 
which supports this picture. In Fig. ||(a) we show the 
low-r region of the PDF from the x = 0.0 and x — 0.25 
samples. Referring to Fig. || we see that these two data- 
sets have relatively narrow Cu-0 bond-length distribu- 
tions. Furthermore, in Fig. ^(a) it is apparent that the 
peak position has shifted due to the change in the average 
Cu-0 bond-length with doping, as expected. The differ- 
ence curve below the data shows that the two data-sets 
are quite different due to the significant structural dif- 
ferences. In Fig. ||(b) we show the intermediate x = 0.1 
data-set plotted as open circles. This peak is centered at 
a position shown by the dashed line which is intermedi- 
ate between the positions of the x = 0.0 and x — 0.25 
data-sets. Referring to Fig. ||we see that the Cu-0 bond- 
length distribution is relatively broad at this composition. 
Plotted on top of the x = 0.1 data-set in Fig. |3|(b) as the 
solid line is the PDF obtained by taking a linear com- 
bination of the x = 0.0 and x = 0.25 data-sets in the 
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FIG. 3. (a) PDFs from the 10 K data from the x = 0.0 
(open circles) and x = 0.25 (solid line) samples. The dif- 
ference is plotted below. The dashed lines indicate expected 
uncertainties due to random errors, (b) PDF from x — 0.1 
data at 10 K (open circles). The solid line shows the PDF 
obtained by making a linear combination in a (1:1) ratio of 
the PDFs shown in (a). See the text for details. 



1:1 ratio, without rescaling the data at all. The differ- 
ence curve is shown below. The good agreement clearly 
demonstrates that the observed PDF peak position and 
broadening of the x = 0.1 data-set is entirely consistent 
with there being an underlying bimodal bond-length dis- 
tribution consistent with heavily doped and undoped re- 
gions of the CuC>2 plane. We can infer from this analysis 
that the difference in the bond-lengths is ~ 0.024 A which 
is the difference between the average bond-lengths of the 
x = 0.0 and x — 0.25 samples. 

Finally, we note that the PDFs from underdoped 
La2- 2; Sr 2 ;Cu04 are consistent with the presence of CuOg 
octahedral tilt disorder in the samples. In an earlier pa- 
per we showed that the measured PDFs could be well 
explained by a local structure which contains a mixture 
of large and small octahedral tilts. We note here that 
we also have evidence from the PDFs for the presence of 



tilt-directional disorder in the sense that there is a mix- 
ture of (100) ("LTO") and (110) ( "LTT" ) symmetry tilts 
present in the local structure. This will be reported in 
detail elsewhere pjj . 

To summarize, we have presented evidence from neu- 
tron diffraction data which strongly supports the idea 
that doped charge in the CuC-2 planes of superconduct- 
ing La2-xSr2;Cu04 for < x < 0.15 and at 10 K is 
inhomogeneous. For doping levels of x — 0.2 and above 
the charge distribution in the Cu-0 plane becomes ho- 
mogeneous. This presumably reflects a crossover towards 
more fermi-liquid-like behavior in the overdoped regime. 
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